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1. Introduction 

Single ionization by electron impact (or (e,2e) collision) at large energy-momentum 
transfer in the nearly Bethe-ridge kinematics, where the recoil ion momentum is 
small compared with the transferred one, constitutes a powerful spectroscopic tool for 
exploring electronic structure of atomic systems [TJ El H] • This (e,2e) method is often 
referred to as the Electron Momentum Spectroscopy (EMS) [HI2|- Theoretical grounds 
of the EMS exploit domination of the lowest-order term, such as given by the plane 
wave Born or impulse approximation (PWBA or PWIA), when treating the ionization 
amplitude perturbatively and using plane waves for description of the incident and 
outgoing electron states. Formally, the ratio of the higher-order and lowest-order Born 

— 1/2 

terms behaves as oc E , where E is the incident electron energy, and hence one might 
expect the lowest-order term to prevail at high values of Eq. However, on the energy 
shell the plane wave higher-order terms are given by divergent integrals j5j, which are 
due to so-called Coulomb singularities of the transition operator [HIE IB]- Therefore, 
for drawing conclusions on the validity of the PWBA model, one should cope with the 
problem of divergences of the corresponding perturbation series. 

The Born series follows from the Lippmann-Schwinger equation which has a 
noncompact kernel. However, in the case of short-range potentials, the noncompactness 
of the kernel does not prevent from calculating matrix elements of the Born series. 
For the purpose of a mathematically correct formulation, Faddeev proposed to 
reduce the Lippmann-Schwinger equation to a system of three coupled equations with 
compact kernels EH|. This system leads to the Born- Faddeev series where two- 
particle amplitudes depending on three arguments (the relative-motion energy and 
incoming/outgoing momenta) appear instead of the corresponding potentials. The 
higher-order Born-Faddeev terms contain the two-particle amplitudes which in general 
are off-shell, i.e. their arguments do not necessarily obey the energy-momentum relations 
for free particles. The lowest-order Born-Faddeev term, given by a half-on-shell two- 
particle amplitude, yields the usual PWIA. 

In the case of Coulomb potentials Faddeev's reduction does not offer mathematical 
advantages, since the kernel remains noncompact [SJEII|- This is a consequence of the 
fact that the plane wave states do not obey the correct asymptotic behavior for the 
Coulomb breakup. Thus, the resultant Born-Faddeev series contains divergent terms 
and, like the PWBA case, the validity of the PWIA model is questionable. At the 
same time, plane waves are usually a convenient and handy mathematical tool for 
calculating the Born and Born-Faddeev series. And as far as the corresponding higher- 
order terms diverge on the energy shell, one must resort to a regularization procedure 
which removes the artificial, unphysical problem of divergences. Further, one must find a 
relation between the regularized perturbation series and the exact ionization amplitude. 
Only then, a perturbative treatment of the ionization amplitude can be developed in a 
physically consistent manner. 

The objective of this work is to revise the traditional PWBA and PWIA models 
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in the light of the fact that these first-order theories ignore the divergence problem 
associated with the higher-order contributions to the ionization amplitude. The 
present theoretical analysis is based on the results of 0, where a practical recipe for 
regularization of the plane wave Born series was proposed, and those of Shablov et 
a>l El IH|, who established a relation between the exact ionization amplitude and 
the unphysical plane wave Lippmann-Schwinger one in the on-shell limit. We proceed 
from the Born series, since, as mentioned above, Faddeev's reduction is not of benefit in 
the case of Coulomb potentials and hence may lead to methodological confusion. After 
regularizing the Born series, one can formulate renormalized analogs of the traditional 
PWBA and PWIA models. These analogs resemble the traditional PWIA model with 
a modified Gamow factor. As shown below, the modified Gamow factor depends on the 
choice of a regularization procedure and therefore it is not uniquely determined. This 
feature permits, in principle, to choose such version of the Gamow factor that efficiently 
incorporates higher-order effects ignored by the traditional first-order models. The above 
points are addressed in the consideration that follows and are illustrated with numerical 
results. 

Section 121 of this paper delivers a general formulation for the (e,2e) reaction on an 
atom. In section El specific Born approximations are formulated using a regularized 
Born series. Section 0] is devoted to the impulse approximation theory in the context 
of the regularization formalism. The numerical results for the case of electron-hydrogen 
ionization are presented and discussed in section El and the conclusions are drawn in 
section El The atomic units (au) e = K = m e = 1 are used throughout unless otherwise 
stated. 

2. General formulation 

We specify the momenta of the incident, scattered and ejected electrons by feo, k s and 
k e , respectively. The corresponding energies are denoted by E , E s and E e . The initial 
atomic and the final ionic states are specified by their respective wavefunctions $f , 
where Z designates the nuclear charge, and energies £», £/. The rate of the (e,2e) 
reaction is characterized by the triple differential cross section (TDCS) 



Here the directions of the outgoing electron momenta are specified by the solid angles 
Q s and Q e . The amplitude T s (T e ) corresponds to the situation where the scattered 
electron has the momentum k s (k e ). The so-called capture amplitude is ignored, 
since we consider such kinematical regimes where the capture of the incident electron 
accompanied with ejection of two atomic electrons having the momenta k s and k e is 
negligible. In equation (JTJ) a sum (average) over unresolved ionic (atomic) states is 
assumed. 

The amplitude is given by (below we focus on the amplitude T s omitting its index) 




(1) 




(2) 
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where Vi is the potential between the incident electron and the atom. The initial 
asymptotic state \k $f) = |fe ) <8> |$f), where |fc ) is the plane wave state for the 
incident electron, satisfies the Schrodinger equation 

(H-V i -E)\k o $f) = 0, 

where H is the full projectile-atom Hamiltonian and E is the total energy: 

E = E s + E e + Sf = Eq + £j. 

The total scattering state \^fj(k s ,k e )) takes account of all interactions between the 
final-state fragments. It satisfies the Schrodinger equation 

(H-E)\Vj(k s ,k e )) = 

and obeys the proper Coulomb asymptotics, whose specificity is due to a long-range 
character of the interactions between the final-state fragments. 

2.1. The plane wave Born series 

To avoid confusion, we will use a tilde for marking the Lippmann-Schwinger analogs of 
the physical quantities that have been introduced in the preceding subsection (such 
as the amplitude T and the total scattering state \tyj(k s ,k e ))). In the context 
of the Lippmann-Schwinger formalism which employs plane wave states for treating 
asymptotically free particles, the total scattering state is sought as a solution to the 
equation 

\*j{k„ fe e )) = Ifc.fce^- 1 ) + G (E)V\^j(k s , fe e )), V = V s + V e + V se , (3) 

where V s ,V e , and V se are the electron- ion and electron-electron potentials, respectively. 
The Green's operator Gq(E) is given by 

G (E) = (E-H + V -iO)- 1 . (4) 

The final asymptotic state \k s k e Q Z f~ l ) = \k s ) ® |fe e ) ® where \k s ) and |fe e ) are, 

respectively, the plane wave states for the scattered and ejected electrons, satisfies the 
Schrodinger equation 

(H- V - E)\k s k e $ z } - 1 ) = 0. 
Equation © can be presented in the equivalent form 

|*7(fc„fc e )> = [l + G-{E)V]\k s k e <$> z f l ) (5) 

where G~~{E) = (E — H — i0) _1 is the full Green's operator. 

Substitution of (JSJ) into (J2J) generates the plane wave perturbation series 

oo 

f = J2T {n) where f ^ = (fe $f {V^Gq (E)V] n \k s k e $ z (6) 

n=0 
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which is traditionally referred to as the Born series. Here the n = term amounts to 
the usual PWBA 

f(p) = t pwba = {k ^z \ V .\ kske <s>z-iy ( 7) 

All other (higher-order) Born terms T^ n ' and the sum of these terms are divergent 
(see jS] and | Appendix A] for details). These observations are not novel and are due to 
a well-known fact that the plane wave states do not obey the asymptotic conditions 
peculiar to break-up Coulomb scattering. Thus, equation © gives an unphysical result 
and hence it can not be employed for a perturbative treatment of the correct, physical 
amplitude (J2J). In addition, the validity of the traditional PWBA model (J7J) can not be 
based upon equation (JHJ). 

2.2. The off-shell Born series 

The situation changes if we consider the Born series © off the energy shell, i.e. when 
E ^ E s + E e + e f . Setting E - E s - E e - e f = A > 0, we have 



oo oo 



f (A) = < n >(A) = T PWBA + J^T (n) ( A ) (T (n) (0) = fW), (8) 



n=0 n=l 



where all the higher-order Born terms are finite p2| (see also |Appendix A ). If the value of 



the off-shell parameter A approaches the on-shell case A = 0, the Born series (JHj) exhibits 
a typical Coulomb singularity T(A) oc A - " 7 El CD], where r\ is a total Sommerfeld 
parameter: 



i — 1> 



n = r] s + rj e + n se (n s = -k s \ r] e = -k e 1 , and r] se = \k s - k e \ 
The physical amplitude (J2J is derived from (jHJ) in the following manner |7| IS} HUj: 

T = eXp( 7^~ iA) lim A-f (A), (9) 
r(l + i?7) a->o K h K ' 

where 

A = r] s \n(2kl) + r7 e ln(2A; e 2 ) + r/ se ln |fe s - fc e | 2 

is the Dollard phase ^T] . Note that in the on-shell limit A — >• the divergent factor A" 7 
compensates for the singularity of T(A). In the next section we show how the calculation 
scheme based on equation @ can be implemented in practice through regularization of 
the Born series. 

3. Born approximations 

For taking the on-shell limit A — > in equation (|9*|) . it is convenient to have A _1J? factored 
out in the off-shell Born series (jHJ. This can be fulfilled by means of a regularization 
procedure (see jAppendix B ). As a result, we obtain the amplitude (jOJ) in the factorized 



form 

oo oo 

t = nf n , f n = J2 t% ] = T PWBA + f n > ( 10 ) 



71=0 n=l 
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where T-£ is the regularized on-shell Born term and 1Z is a regularization functionf, 
such that 1Z = 1 if rj s = r\ e = r\ se = 0. 

In contrast to ©, the result (fTUjl allows to develop a perturbative treatment of 
the physical amplitude. However, when truncating the regularized on-shell Born series 
T-ji to a finite number of terms, there is an uncertainty associated with the choice of a 
regularization procedure. Namely, while 

nf n = K'f n ,, 

where VJ and T^i are due to an alternative regularization procedure (see Appendix B 
for details), in general 

N N 
n=0 n=0 

In particular, the lowest-order approximation (N = 0) to the amplitude (fTU|) assumes 
the form 



T PWBA = n f(0) = nT PWBA (n) 

that depends on the choice of a regularization procedure/function. These observations 
are reminiscent of the situation that one encounters in quantum electrodynamics, 
when regularizing a series of Feynman's diagrams. Drawing an analogy with that 
situation, one can speak of a renormalization group which is formed by the regularization 
procedures in the present case. Therefore we will refer to equation (fTTjl as a renormalized 
PWBA (RPWBA). 

In this work we develop the Born series for the physical amplitude in a manner 
similar to equation (0), that is 

oo 

T = J2 T(n \ ( 12 ) 

n=0 

where is given by the nth-order term of the Maclaurin series expansion of the exact 
amplitude T with respect to the two-particle Sommerfeld parameters rj s , rj e , and rj se . It 
can be deduced from (fjT)j) that 

T (n) = fW^ ( 13 ) 

where T^li corresponds to the specific regularization procedure which yields TZ = 1 for 
any values of rj s , rj e , and r/ se (see |Appendix B ). 

According to ([EH) , the lowest-order term of (|T2|) amounts to the usual PWBA 
result ©: 

T (0) = t pwba = ( ko $f\ V .\ kske <f>z-ly (14) 

The plane wave second Born approximation (PWB2), which is usually of practical value 
for estimating the applicability of the PWBA, is then given by 

yPWB2 _ rp(0) _|_ rp(l) _ /jiPWBA _j_ ji(l) 
] Hereafter 1Z refers to an arbitrary regularization procedure/function unless otherwise specified. 
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Using (JHJ), this is equivalent to 

t pwb2 = ^ _ ™7 _ iA+ ^ t pwba + hm[T«(A) + vrjT PWBA In A], (16) 
where 7 = 0.577216 is the Euler constant. 

4. Impulse approximations 

4-1. The usual theory 

Near the Bethe ridge the (e,2e) process can be modelled as a binary encounter between 
the projectile electron and the electron that is ejected from the atom pQ. The 
PWIA theory formulates this picture mathematically. Within the Lippmann-Schwinger 
approach (jHJ) it is expressed as follows: 

00 

yPWIA = where f{n) = (fe $f I V se [Gq (E) V se } n \ kske^j' 1 ) . (17) 

n=0 

As in the case of ©, the n > 1 terms and the sum of these terms are divergent. We 
can remedy this defect in order to obtain a correct, convergent counterpart of (|17|). 
For this purpose we use the result (fTUj). assuming that Vi = V = V se and, accordingly, 
Vs = Ve = 0. Thus one derives the physical counterpart of (fT7|) as 

t pwia = n j nF . f ^ = TF . f ( q ) (7 ^ e = n\ Vs=Ve=0 ), (is) 

where q = k s + k e — fc is opposite to the recoil ion momentum and Fif(q) = ($f |$J _1 qr) 
is the so-called structure amplitude lj. t-ji is the regularized Lippmann-Schwinger half- 
on-shell amplitude for ee-scattering [3] and r is the exact half-on-shell ee-scattering 
amplitude [P^ITS]: 

r = r{k - \{k s + fc e ), \{k s - fc e ); [\{k s - k e )} 2 }. 

The result (fTHj) is irrespective of the choice of a regularization procedure, since TZ se t-ji = 
7Z' se tft'. Note that equation (fTHjl is a corner-stone of the usual PWIA theory of (e,2e) 
reactions on atoms PUI2]- 

4-2. A renormalized theory 

While there is nothing wrong in the above derivation of the traditional PWIA result, it 
can be noticed that at the starting point (fTTj) we have focused on the purely electron- 
electron part of (JBI) and thereby we have ignored the remaining part which is also 
divergent. Methodologically, one should apply the binary-encounter approximation to 
the regularized Lippmann-Schwinger amplitude T-^ (see (jlUj) ) which, in contrast to 
is free of divergences. The electron-electron component of Tn is given by the regularized 
Lippmann-Schwinger amplitude in for ee-scattering (see (fTSjO . Using (fTU|) and (JTHJ), the 
PWIA to the physical amplitude T is then given by 

T£ WIA = Ki n F lf {q) = (n/n se )rF if (q). (19) 
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Since in general 1Z/1Z se ^ TZ'/TZ' se , the PWIA amplitude (fTU|) is determined by the choice 
of a regularization procedure. In particular, the specific case 1Z/1Z se = 1 amounts to the 
traditional PWIA (j!8j) . As far as the choice of a regularization procedure is a matter of 
taste, equation ()19|) offers an infinite number of alternatives. This fact has the following 
general consequence: there is no a priori PWIA model in the case of (e,2e) reactions on 
atoms. By analogy with the RPWBA case (jllj) . the approximation (fK?|) will be referred 
to as a renormalized PWIA (RPWIA). 

4-3. TDCS and the Gamow factor 

Within the traditional PWIA model (JIBJ) the TDCS (0) is expressed as P 
d 3 a PWIA k,k e g(n se ) 

k - k. 



\k 


k s \ 4 






\k 


~~ k s | 4 




k s | 


\k 


ke\^ 


\ko 


h |2 
"'el 



1 + 71 m _ 71 n? cos Vse in 



E( av ) o 
\Mq)\ 2 , (20) 

where denotes the average over initial-state and sum over final-state degeneracies, 

and Qirjse) is the so-called Gamow factor [TJ^ 



g{n se ) = | exp(-i7rr7 se )r(l - ir] se )\ 2 = J^**L \ ' ^ 

In the case of the RPWIA model the expression for TDCS depends on the explicit 
form of the regularization function TZ. In this work we inspect the following form: 

1Z = exp(— ^Trr])T(l — ir/) and consequently lZ se = exp(— ^Trrj se )T(l — ir] se ). (22) 

Substitution of (HHJ) and (J22J) into © yields 
dV RPWIA k s k e g( v ) 

\ko — k s 



\k 


~ k s |^ 






\k 


~ k s | 4 


\ko 


k s \ 2 


\k 


k e \^ 


|feo 


h |2 
"'el 



1 + ti m ~ 71 n? cos Vse In 



k - k 



x 2. 



\F if (q)\ 2 , (23) 



with the modified Gamow factor 

«•»> - ^TT' ' 24 » 

As can be deduced, the RPWIA result (|23|) differs from the traditional PWIA 
one ()20|) only in the expression for the Gamow factor. In contrast to (}2~T]) . the modified 
Gamow factor (|24|) treats the final-state particle pairs on equal footing. It should be 
noted that the same result as P^j) is obtained using the first-order model introduced 
by Shablov et al [7] on the basis of the formalism of regularization Coulomb operators, 
and it can be traced down to the idea of effective charges proposed by Peterkop [15j. 
The influences of the traditional and modified Gamow factors on TDCS are examined 
in the next section, where we present the corresponding numerical results. 
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Using the RPWBA model and equation (|22|) . in the binary-encounter 

approximation (Vi = V se ) we get 



d 3 



a 



RPWBA 



d 3 



a 



PWBA 



dQ s dQ e dE e y " dtt s dn e &E e 
where the traditional PWBA result is given by pQ 

d 3 (T PWBA Kk P 1 

1 + 



(25) 







\k - 




\k 


— k e \^ 


|feo - 


fc e | 2 



E 



(av) 



dn s dVl e dE e 2ir 3 k \k - fc s | 4 

It can be noticed that in symmetric kinematics (i£ s = E e , \k — k s 
is usually the case of the EMS experiments (see 13 H] and references therein), the 
TDCS given by equation (|2"K|) is identical to that given by equation (J22J). 



^/(<?)| 2 - (26) 
= \k — k e \), which 



5. Numerical realization 



In this section we present numerical results for an archetypical case, namely the (e,2e) 
reaction on a hydrogen atom. The symmetric setup E s = E e and 9 S = 9 e = 45° is 
inspected, where the polar electron angles 9 S and 6 e are measured with respect to the 
direction of the incident electron momentum. In the considered setup the TDCS is 
usually studied as a function of q = \q\ (see (fTTjl) whose value is varied in noncoplanar 
geometry by varying the value of the relative azimuthal angle A<ft se = cf) s — <f) e , where <p s 
(4> e ) is the azimuthal angle of the scattered (ejected) electron. The minimal value of q 
corresponds then to the symmetric coplanar case (|A0 se | = tt), where the incident and 
outgoing electron momenta are in the same plane. 

We focus on high energy values (E s — E e > 1 keV) which closely meet 
the Bethe-ridge and binary-encounter criteria. In the absence of the corresponding 
EMS measurements, particularly those performed on an absolute scale, the results 
of the traditional and renormalized first-order treatments are compared with those 
of the PWB2 calculations and those of the Brauner-Briggs-Klar (BBK) model ^E], a 
representative of nonperturbative treatments. The PWB2 calculations have been carried 
out in accordance with ([16)1 . where the on-shell limit has been taken analytically using 
the regularization procedure described in |Appendix B| The regularized PWB2 integrals 
have been performed numerically following the method developed in [T7] . 

In the BBK model the final state in (j2J) is given by 

|*B BK (fc s ,fc e )) = \<p-(k,)) <g> \<p~(k e )) ® \x~(k se )) (k se = |(fe s -fc e )), (27) 

where \(p~(k s )) and \ip~(k e )) are the Coulomb waves describing outgoing electrons 
moving in the field of the proton. The electron-electron correlation factor \x~(k se )) 
is determined through 

\k a ) ® \k e ) ® \x~(kse)) = \4>~(k s ,k e )), 

where \ip~(k s ,k e )) is the exact scattering state in the absence of the electron-proton 
interactions. The wave function (|2"7j) has the correct asymptotic behavior and gives 
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correct results in the limiting situations, where (i) the charge of one of the final-state 
particles is switched off and (ii) the electron-electron interaction is absent. Note that 
the BBK model yields the regularization function ()B.7)) as 

R BBK = exp(-±7rr/)r(l - ir} s )T{l - in e )T(l - in se ). (28) 

The corresponding Gamow factor is then given by 

2tCT} s 27T77 e 27T1] se 

Geek = QWWiVse) = _ 1 ^ _ 1 e2m?ae _ r ( 29 ) 

Like (|24j) . this Gamow factor also treats the final-state particle pairs on equal footing. 

The numerical calculations using the traditional PWBA and PWIA models 
have been performed in accordance with equations (|26|) and (|20|) . respectively. As 
remarked in subsection 14 M\ in the case of symmetric kinematics the RPWBA (J25|) and 
RPWIA ()23|) models are equivalent and therefore below the corresponding numerical 
results referred to as RPWBA/RPWIA. For a hydrogen target, the structure factor 
entering equations (JUJ), ((2HJ), (1231). and lEfil) is 

E (av V,/(9)r = i^(^r, (so) 

where <pi s {q) is the Is state momentum-space wave function. 



5.1. Numerical results and discussion 

Figure Q shows the numerical results for the symmetric noncoplanar kinematics 
utilized in the recent (e,2e) measurements on helium j^j. It can be seen that the 
RPWBA/RPWIA results are substantially larger in magnitude than the PWIA ones, 
which are the smallest in magnitude. This feature indicates an appreciable role of 
the choice of the Gamow factor in the kinematics under consideration. Interestingly, 
the PWBA and BBK results are close to each other both in magnitude and in shape. 
Though the BBK model is not exact, it takes into account those higher-order effects 
that are entirely neglected by the PWBA treatment. One might thus conclude that the 
higher-order contributions to the TDCS are subsidiary in the present case. However, 
this conclusion needs experimental verification, since marked discrepancies between the 
PWBA and PWB2 results are observed in figure H It should be remarked that the 
developed renormalized first-order theories give an opportunity to fit the exact TDCS 
by the proper choice of the Gamow factor. For example, setting (cf equation (|24j) ) 

G(v) = 1 

in equations (f2*3|) and (J2*5|) . one obtains the traditional PWBA result and thus, as can 
be seen in figure ^ rather well reproduces the BBK results for the present kinematics. 
And the PWB2 results are satisfactorily reproduced in magnitude using the modified 
Gamow factor (J24"j) . 

In figure H the traditional PWBA and PWIA and the RPWBA/RPWIA values 
are presented on a logarithmic scale. It can be seen that all three models exhibit 




Figure 1. TDCS as a function of the absolute value of the recoil ion momentum in 
the symmetric noncoplanar kinematics (E s = E e = 1000 cV, 9 S = 9 e = 45°) of the 
recent (e,2e) experiments Un- 



practically identical shapes for the TDCS. This feature is due to the almost constancy 
of the Sommerfeld parameters r] se and 77, and hence of the Gamow factors (|21|) and ([24)1 . 
in the involved kinematical region. It should be noted that the value of the Gamow 
factor ()29|) is nearly the same as that of (J2"4"j) . This observation can be explained by 
the feature that the first-order expansions of (j2IJ) and (j2HJ) with respect to r] s , r] e , and 
r] se are identical. Another important observation is that, in contrast to the traditional 
Gamow factor (j2~T| . the modified Gamow factor (|24J) enhances the magnitude of the 
TDCS with respect to that in the conventional PWBA case. At the same time, all 
three first-order models are practically equivalent for description of the corresponding 
symmetric noncoplanar measurements provided that the latter are performed on a 
relative scale (see, for instance, JHI)- This conclusion follows from the fact that in the 
case of the inspected first-order models the shape of the TDCS (the so-called momentum 
profile [H El El ) is almost fully determined by the structure factor (JHU)) . 

The effect of the Gamow factor on the magnitude of TDCS in the considered 
geometry is shown in figureEl As can be deduced from the figure, the RPWBA/RPWIA 
and the traditional PWIA values for the TDCS rather slowly converge to each other 
and ultimately to the traditional PWBA value upon the increase of the incident electron 
energy E Q . For E ~ 10 keV, the RPWBA/RPWIA results are larger in magnitude than 
the traditional PWIA ones by a factor of ~ 1.5, and even in the region Eq ~ 100 keV, 
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Figure 2. The first-order results for the TDCS. The kinematics is the same as in 
figure CI 

where one might expect relativistic effects to come on the scene, the relative difference 
between the results in magnitude is about 10%. These findings are in discord with an 
intuitive physical picture which assumes practical equivalence of the first-order models 
in the kinematics that approaches the classical ridge of a billiard-ball collision 2j. The 
discrepancy can be explained by the following factors: (i) the long-range Coulomb forces 
between the colliding electrons, as opposed to the contact-like forces between classical 
billiard balls, and (ii) the presence of the Coulomb field of the ion (this factor is relevant 
only to the RPWBA/RPWIA model). 

6. Summary and conclusions 

In summary, we have considered the electron-atom ionization process at large energy- 
momentum transfer and near the Bethe ridge. Proceeding from the plane wave 
Lippmann-Schwinger amplitude, which diverges on the energy shell, we have regularized 
the corresponding Born series. On this basis we have developed the correct, physical 
Born treatment whose lowest-order term amounts to the conventional PWBA. The 
RPWBA model has been formulated, which depends on the employed regularization 
procedure. We have shown that the PWIA can not be uniquely determined and therefore 
we have introduced the RPWIA model, in which the Gamow factor is determined 
by the choice of a regularization procedure. The numerical results for the symmetric 
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1 10 100 
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Figure 3. The traditional Ij21|) and modified (|24|l Gamow factors as functions of the 
incident energy E in the symmetric coplanar case (E s — E e , 9 S = 9 e = 45°, and 

\A(j) se \ = 71"). 

noncoplanar kinematics have been presented. It has been demonstrated that even at 
high electron energies (~ 10 keV) the magnitude of the TDCS is very sensitive to the 
choice of the Gamow factor. 

We might expect the results of this work to be primarily important for theoretical 
treatments of absolute (e,2e) measurements in the nearly Bethe-ridge kinematics at large 
energy-momentum transfer and for further development of the EMS method pQ. The 
present theoretical consideration can be generalized to the case of ionization of an atom 
by a charged-particle impact, for example, to the cases of positron- and proton- atom 
ionization. Using the formulated RPWBA and RPWIA models with a properly modified 
Gamow factor, one can efficiently take into account the higher-order effects ignored 
by the traditional PWBA and PWIA models. In this connection, it should be noted 
that a consistent, rigorous treatment of the higher-order contributions to the ionization 
amplitude is realized by means of the developed plane wave Born series ()12)) . which 
has an apparent advantageous feature: the value of any of its terms is, by definition, 
irrespective of the choice of a regularization procedure. 
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Appendix A. Divergencies of the Born series 

To elucidate the origin of divergencies, let us examine the electron-electron part of the 
n = 1 term in © 

se ~ J (2n)3(2Ti)3E-p2/2-pl/2-e f -i0 [PsPel sel 3 e) 
dp {k^f\Vi\k s -p,k e + p^ z f - l )ATi 

(A.l) 



(2tt) 3 (k s -k e )-p-p 2 -i0 p 2 ' 
It can be seen that the integrand has a pronounced singularity at the point p = which 
physically corresponds to elastic rescattering in the forward direction. This feature 
makes the integral (|A.1J) divergent. The same argument applies to the electron- ion 
parts of the n = 1 term, since for the static matrix elements of V s and V e we have: 

4:71 

(k s ± P ,$ z f - 1 \v s \k s $ z f - 1 ) = {k e ±p,<$> z f - l \V e \k e <$> z f - 1 ) ~ -— (p -►()). 

Clearly, in the case n > 2 we encounter even stronger divergencies because the 
corresponding Born terms contain multiple elastic rescattering in the forward direction. 
In the case of (|S]L the electron-electron part of the n — 1 term takes the form 

(cf flu)) 

fW(A) f dp Wl^ fc '-^ + P-y) ^ 
seK 1 J (2vr) 3 A + (k 3 - k e ) ■ p - p 2 - iO p 2 ' 1 ' ' 

where the singularity in the integrand at the point p = is removed (note that 
dp = p 2 dpdQ p ) and thereby the integral does not diverge. Specifically, in the case 
A — ► we have [Sj 

f«(A)~-ir /se T PWBA lnA. 
Appendix B. Regularization procedure 

Below we describe a possible recipe for factoring out A _ir? in equation (JHJ) and then taking 
the on-shell limit in equation It consists in presenting the Green's operator (j3J) in 
the form 

G o (*0 - 2. J J^J^E-p 2 /2-p 2 e /2-e f/ -iO ~ Q ° {E) + ^ {E) ' (B ' 1} 
where Qq(E) is the regularized Green's operator and 

\Z-1\ 1 kZ-1\ 



^ [ > ZjfJ (2vr)3 (2vr)3 E - p 2 /2 - p 2 /2 -e f ,-iO 



(B.2) 
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is the Green's operator component which is responsible for divergencies. Using (|B.1|) 
and (jB.2|) . we obtain the off-shell Born series (JBJ) in the factorized form 

T(A) = P R (A)f R (A), (B.3) 



where 



with 



f R {A) = £f< ft) (A) = T PWBA + J2t ( r\A), 

n=0 n=l 

CO oo 

Ph(A)=53pW(A) = 1 + ^pW(A) > 

n=0 n=l 

i** {A) = (fco*f |y i K(B)V]"|fc,fc B *f- 1 >, 
i#'(A) = E / ^^<fw.«f- , |[GJ(B)VT|M.«f- 1 >. 



Taking into account that 



l + J2[Go(E)V\ n = l + G-(E)V 



n=l 



and using ©, we get 

Um P H (A) = J]/ ^^(p^-^j^,*.)). (B.4) 

The on-shell limit f)B.4|) does not exist, since the Lippmann-Schwinger total scattering 
state \^fj(k s , k e )) is not physical. Using the results of Shablov et al 0|B], we deduce 
that 

P R (A 0) = A^ exp(ivrr/ + L4)T(1 + if?) 

x E / ^^^r\^k e )). (B.5) 
Inserting ()B.3|) into (jOJ) and using (|B.5|) . we obtain the physical amplitude as 



T = PT fi , where f R = ^ T^ n) = T PWBA + £ t£°. (B.6) 

n=0 n=l 

Here = liniA->oPjj^(A) is the regularized on-shell Born term. The regularization 
function R is given by 

R = ^(5(r s )5(r e )$fr 1 |*7(fc s , fe e )), (B.7) 

/' 

where 5(r) designates Dirac's delta function. The function ()B.7|) has the following 
obvious property: R = 1 if f] = 0, where i) = (r) s , r] e , rj se ). 

Note that the above recipe for regularization is only one among an infinite number 
of possible regularization procedures and, in general, one obtains different regularization 



Electron-atom ionization near the Bethe ridge 



16 



functions for different regularization procedures. It means that while the product 
RT R = T (jB.6|) is an algoristic quantity, the factors R and T R are not, i.e. RT R = IZTn, 
where TZ {TZ = 1 if f] = 0) and T-r corresponds to an alternative regularization procedure. 
For example, there is such regularization procedure that yields TZ — 1 for any value of 
T). To illustrate this statement, we expand R in the Taylor series with respect to the 
components of fj and notice that T R cx fj n . We have 



R Z_> n j fl " Z^Z_^ (n-m)! fi 

n=0 m=0 n=0 m=0 



Y,fH = f n=u (B.8) 



where it is supposed that the operator 



T] ■ V« =0 = % 7^— + ?7e 7^— + Vs 



9VsJfj =0 \dVeJfj= \dVse /n=0 

acts only on the regularization function R, and (i) ■ V^=o)°-R = -R|^=o = 1- As can be 
deduced, oc r and f^ =1 = f R 0) = T PWBA . 
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